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ABSTRACT 

Hubble Space Telescope observations of distant clusters have suggested a steep increase 
in the proportion of SO galaxies between clusters at high redshifts and similar systems 
the present-day. It has been proposed that this increase results from the transformation 
of the morphologies of accreted field galaxies from spirals to SOs. We have simulated 
the evolution of the morphological mix in clusters based on a simple phenomenological 
model where the clusters accrete a mix of galaxies from the surrounding field, the spiral 
galaxies are transformed to SOs (through an unspecified process) and are added to the 
existing cluster population. We find that in order to reproduce the apparently rapid 
increase in the ratio of SO galaxies to ellipticals in the clusters, our model requires 
that: (1) galaxy accretion rate has to be high (typically, more than half of the present- 
day cluster population must have been accreted since z ~ 0.5), and (2) most of the 
accreted spirals, with morphological types as late as Scdm, must have transformed to 
SO's. Although the latter requirement may be difficult to meet, it is possible that such 
bulge-weak spirals have already been 'pre-processed' into the bulge-strong galaxies 
prior to entering the cluster core and are eventually transformed into SO's in the cluster 
environment. On the basis of the evolution of the general morphological mix in clusters 
our model suggests that the process responsible for the morphological transformation 
takes a relatively long time (~ 1-3 Gyr) after the galaxy has entered the cluster 
environment. 
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1 INTRODUCTION 

The morphologies of galaxies remain one of the key observa- 
tional characteristics used to classify and differentiate classes 
of galaxies. It is known from observations of the local Uni- 
verse that galaxy morphology depends on environment, with 
the galaxy population within clusters dominated by early- 
types and spirals being found predominately in lower-density 
environments. Even within the clusters, the morphologi- 
cal mix changes with local galaxy density (and hence with 
cluster-centric distance) - producing a morphology-density 
relation - as shown by Dressier (1980, D80). Moreover, sim- 
ilar segregation behaviour has been seen in Hubble Space 
Telescope (HST) imaging of ten clusters at higher redshifts, 
z = 0.37-0.56 by Dressier et al. (1997, D97). These results 
confirm galaxy morphologies as one of the most sensitive 
tracers of environmental influences currently known. 

One of the most interesting results to come out of the 
morphological studies of distant clusters has been the claim 
by Dressier et al. (1997) that the SO-to-elliptical (S0/E) ra- 
tio is a strongly decreasing function of redshift: by a factor 



of ~ 5 between z ~ and 0.5 (although see Andreon 1998; 
Fabricant, Franx & van Dokkum 2000). This contrasts with 
the increasing number of blue, spiral galaxies seen in these 
regions at higher redshifts; the origin of the Butcher-Oemler 
effect (Butcher & Oemler 1984; Dressier et al. 1994; Couch 
et al. 1994, 1998; D97). D97 therefore proposed that most 
of the cluster SO's are formed relatively recently by the mor- 
phological transformation of spirals (see also Poggianti et al. 
1999; Fabricant et al. 2000). These spiral galaxies are con- 
tinuously supplied by accretion from the surrounding field in 
the course of the assembly of the cluster. Couch et al. (1998, 
C98) and Fasano et al. (2000) have added estimates of the 
S0/E ratios for three clusters at z = 0.31 and nine clusters 
at 0.09 < z < 0.26, respectively, to bridge the gap in the 
D97 analysis between z = 0.37 and the present-day. These 
studies confirm the trend of a decline in the cluster SO pop- 
ulation from z — out to z £ 0.5. Proposed mechanisms for 
transforming the morphologies of bulge-strong spirals to SO 
galaxies include dynamical interactions, such as harassment 
(Moore et al. 1996, 1999a), while gas dynamical processes, 
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e.g. ram- pressure stripping (Abadi, Moore & Bower 1999), 
will also aid in truncating star formation and hence fading 
the disk component of any spiral galaxy within the cluster 
environment. 

The luminous galaxy populations of rich clusters at z ~ 
are dominated by passive galaxies with early-type mor- 
phologies and apparently old stellar populations, as shown 
by their tight colour-magnitude relation (e.g. Bower, Lucey 
& Ellis 1992; Terlevich et al. 1999). The transformation of 
star-forming spiral galaxy to passive cluster member must 
occur, as the blue, spiral galaxies we see in clusters at z ~ 0.5 
cannot escape and hence their descendents have to reside in 
these regions at the present-day (Bower, Kodama & Ter- 
levich 1998; Smail et al. 1998; van Dokkum et al. 1998). 
Following the work by Smail et al. (1998) and Poggianti 
et al. (1999) on the relationship between the blue galax- 
ies in distant clusters and nearby cluster SO's, Kodama & 
Bower (2001) quantified the analysis by sketching the flow of 
galaxies across the colour-magnitude diagram down to the 
present-day, taking into account the galaxy accretion from 
the surrounding fields. 

Motivated by this success in connecting the photomet- 
ric evolution of the galaxy populations in distant clusters to 
the present-day counterparts within the framework of the 
assembly of rich clusters, we now try to seek to explain the 
evolution of the morphological mix in clusters, and espe- 
cially the decline in the SO population, using a simple model 
for the growth of clusters and the evolution of their galaxy 
populations. 

Since the proportion of spiral to elliptical galaxies in the 
field is much higher than that in clusters at any time from 
z ~ 0.5 down to the present-day (e.g. Driver et al. 1998; 
van den Bergh et al. 2000), continuous accretion of galax- 
ies from the field should be an effective process in chang- 
ing the morphological mix in clusters. If field spiral galax- 
ies are turned into SO's upon accretion on to clusters, then 
the accretion process effectively increases the fraction of SO 
galaxies with time, and might explain the rapid increase 
of S0/E ratio presented in D97. The aim of this Paper is 
to test this hypothesis quantitatively and thereby constrain 
both the galaxy accretion rate and the timescale of the mor- 
phological transformation. Unless otherwise stated, we use 
Ho = 50kms _1 Mpc -1 and q = 0.1. 



2 THE MORPHOLOGICAL MIX IN DISTANT 
CLUSTERS AND THE FIELD 

2.1 The Morphological Mix in Clusters 

To investigate the morphological evolution of the galaxy 
population in clusters we use observations of galaxy mor- 
phologies in three samples of rich clusters spanning z = 0- 
0.5. These come from the survey of local clusters by Dressier 
(1980), where we use the same ten high concentration clus- 
ters as in D97, and two different HST-based studies: Couch 
et al. (1998)'s work on three clusters at z = 0.31 and the 
morphological catalogue of ten rich clusters at z = 0.37- 
0.56 from Smail et al. (1997b), which is the basis of the D97 
analysis. 

To reliably compare the galaxy samples at different red- 
shifts we must apply the same selection criteria to each sur- 
vey. Specifically, when attempting to trace the evolution of 



galaxies between high redshifts and the present day, we have 
to take into account the various degrees of fading they will 
suffer after their star formation ceases, which depends upon 
their star formation histories. For example, the most vig- 
orously star forming late-type spirals in the distant cluster 
samples may fade substantially (<; 1 mag) after they cease 
star formation (see Kodama & Bower 2001) and therefore 
only the brighter examples (My g —21) will appear in a 
local cluster sample limited at My ~ —20. For the simple 
model discussed here we have chosen to ignore the added 
complication of a star-formation dependent magnitude cut 
and simply adopt a fixed magnitude limit for all morpholog- 
ical types. The D80 sample is limited at an absolute mag- 
nitude equivalent to Mv ~ —20.4 (D97). We have therefore 
also chosen to limit the more distant samples from D97 and 
C98 at an absolute magnitude of Mv = — 20 (we note that 
the magnitude limit quoted in D97 for their analysis is incor- 
rect and actually corresponds to Mv = —19). By adopting a 
fixed absolute magnitude cut we will tend to over-estimate 
the proportion of descendents of late-type spirals predicted 
in the local clusters. However, we note that the fraction of 
late spirals (Scdm) is small (7 per cent) and therefore our 
results do not significantly change if we use, for example, a 
My = —21 cut for this population. 

In addition, because of the morphology-density relation 
(D97), it is important to use equivalent areas in the nearby 
and distant clusters when comparing the galaxy populations. 
The morphological mix in the high redshift clusters in D97 
and C98 is calculated within the HST fields, which roughly 
correspond to the inner 130" or 1.0 Mpc (z ~ 0.5) diameter 
regions. For the nearby sample from D80 therefore we use 
the same physical radius to define the morphological sample. 

We show in Figs. 1 and 2 the observed S0/E and Sab/SO 
ratios, respectively, from D80, C98 and D97. For the local 
sample of D80 we plot both the individual clusters as well as 
the average and the scatter for the whole sample. We note 
that observations of all of the ten clusters are consistent with 
the mean value given the large uncertainties of the individual 
cluster measurement arising from sampling errors. 

2.2 The Morphological Mix in the Field 

To successfully model the evolution of the morphological 
mix in clusters we need to know the distribution of mor- 
phologies for the galaxy population they accrete from the 
field. Ideally the morphological mix for this field sample 
would be based on a large, volume-limited sample culled 
from a morphologically-classified faint field galaxy redshift 
survey. Unfortunately, there are no large redshift surveys 
with adequate morphological information and so we have 
instead chosen to use a magnitude-limited morphological 
survey at a depth which should place the majority of the 
galaxies within the redshift range we require. We therefore 
use the morphological classifications from the HST-based 
Medium Deep Survey (MDS, Griffiths et al. 1994) which 
reaches I = 22 and will have a median redshift of z ~ 0.5 
with 75 per cent of the galaxies at z < 0.65, based upon 
the redshift distribution at this depth in Lilly et al. (1995). 
This means that the majority of the sample represent the 
population which will be accreted onto the clusters in the 
redshift range covered by our analysis. The morphological 
mix of field galaxies for this sample is then: E, 0.10 ± 0.01; 
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SO, 0.10 ± 0.01; Sab, 0.25 ± 0.01; Scdm, 0.30 ± 0.01; and 
Irr, 0.25 ± 0.01 (as used by Smail et al. 1997b). The errors 
are based on Poisson statistics. This distribution is consis- 
tent with the smaller sample of van den Bergh et al. (2000) 
in the Hubble Deep Field with spectroscopic redshifts which 
place them at z = 0.25-0.6, as well as the small number of 
field galaxies with morphological and spectroscopic informa- 
tion from Dressier et al. (1999). This then provides us with 
a morphological distribution for the field galaxies which are 
accreted onto clusters at z < 0.5. 

In a similar manner to that discussed above, when com- 
paring the cluster populations at different epochs, we also 
need to consider a variable magnitude limit when estimating 
the properties of the field galaxies which are accreted onto 
the clusters, depending upon their previous star formation 
histories (or crudely the morphological type). However, the 
lack of any large morphologically-classified faint galaxy red- 
shift surveys also means that we do not have the information 
available to reliably apply a variable apparent magnitude cut 
to the different morphological subsamples in the field pop- 
ulation. To quantify the likely error which will arise from 
adopting a single magnitude cut across all morphological 
types we turn to the type-dependent luminosity functions 
available for local galaxy samples* from the 2dF Galaxy 
Redshift Survey (Folkes et al. 1999). Correcting the B-band 
luminosity functions to U-band using for type-dependent 
(B - V) colours from RC3 (Buta et al. 1994), and adopt- 
ing limits of M v = -20 for field E/SOs and M v = -20 
or My — —21 for spirals (Sab+Scdm) we can estimate the 
likely variation in the fractions of E/S0, Sab and Scdm. Us- 
ing a limit of Mv = —21 for the spiral galaxies, we es- 
timate relative proportions of 1:0.90:0.95 for E:Sab:Scdm, 
compared to 1:1.25:1.5 for the Mv = —20 cut (we note that 
this latter ratio is very similar to the value derived above for 
the distant field, using Mv = —20). We conclude therefore 
that by adopting a fixed magnitude limit across all morpho- 
logical types in the field we increase the proportion of spiral 
galaxies in the field by factor of ~ 1.4-1.5 for both early- 
and late- type systems. We will discuss the consequence of 
this effect in 54. 



3 MODELLING THE EVOLUTION OF 

GALAXY MORPHOLOGY IN CLUSTERS 

Using the morphological distributions for the cluster and 
field populations derived in the previous section we devise 
a series of models for the evolution of the morphological 
mix in clusters from z ~ 0.5 to z = 0. These models are 
by necessity simple and ignore many of the details of the 
dynamical processes which may be involved in the morpho- 
logical transformation, such as galaxy harassment (Moore et 
al. 1996; 1999a) and/or ram-pressure stripping (Abadi et al. 
1999). However, the relative importance of these processes 
is still a matter of debate (e.g. Poggianti et al. 1999) and the 
details of the processes are highly uncertain. We therefore 

* While there is no suitable high redshift sample to use, we note 
that there is no evidence for differences in the evolution of the 
bright end of the luminosity functions for colour-selected (red 
and blue) galaxy populations out to z ~ 0.5, at least down to our 
limit of M v = -20 (Lilly et al. 1995). 



Table 1. Summary of the Models 



Models 


A 1 " 

recent 


Morphological 


Ttrans 






Transformation 


(Gyr) 


(a) hi-all-2 


high 


Sab+Scdm -> SO 


2 


(b) lo-all-2 


low 


Sab+Scdm -> SO 


2 


(c) hi-early-2 


high 


Sab ->■ SO 


2 


(d) hi-all-3 


high 


Sab+Scdm -> SO 


3 


(e) hi-all-1 


high 


Sab+Scdm -> SO 


1 


(f) shi-all-1 


super-high 


Sab+Scdm -> SO 


1 



Note — t The three different accretion rate models have 
cumulative galaxy accretion fractions at z = 0.5 and 0.3, 

(A reC ent(2 = 0.5), A r 

cccnt 

(z = 0.3)), of: low, (0.35, 0.2); 
high, (0.5, 0.3); and super-high, (0.7, 0.4). 



adopt a phenomenological approach to test the basic foun- 
dations of the morphological transformation hypothesis for 
the evolution of the cluster SO population. 

We compare the predictions of the various models for 
the evolution of the S0/E and Sab/SO ratios within the clus- 
ters with the observations from D80, C98 and D97. As we 
show below, these comparisons constrain both how large a 
fraction of the cluster population must be accreted since 
z ~ 0.5 and how long the morphological transformation 
takes. 

There are three key parameters in the models: (1) the 
total galaxy accretion rate onto the cluster; (2) the range of 
spiral types (e.g. Sab, Scdm) which transform into cluster 
SO's, and (3) the timescale of this morphological transfor- 
mation. We note that the accretion rate and morphological 
range have almost the same effect on the morphological evo- 
lution in the clusters, since both determine how many of the 
accreted field galaxies are turned into SO's. In contrast the 
transformation timescale controls how long the accreted spi- 
rals are classified as having 'spiral' morphology after enter- 
ing the cluster. The parameters for the models we consider 
are summarised in Table 1. 

The cumulative galaxy accretion rate onto a cluster, 
A cuul (t), is normalised to unity at z = (equivalent to t = 
tuniv, the age of the present-day Universe — 16.6 Gyr in our 
cosmology) as: 

^eum(tuniv) = / A(t)dt = 1, (1) 

Jo 

where A(t) is the normalised accretion rate at a cosmic 
time t. To quantify the amount of recent accretion we de- 
fine A rC ccnt(i), the integrated accretion rate from time t (or 
redshift z) to the present-day (t = i U niv, or z — 0): 

A r ecent (^) — ^cum (^univ) ^cum(^)- (2) 

The actual form for the accretion rate is taken from the 
Extended Press- Schechter theory (Bower 1991), assuming a 
power spectrum with n = —1.5. Based upon the weak lens- 
ing estimates of the masses of the clusters used in D97's anal- 
ysis (Smail et al. 1997a) we adopt a mass for the present-day 
cluster of ~ 10 15 M Q , or 100 M* . Since the galaxy accretion 
history is strongly dependent on the threshold mass for the 
substructure within the clusters (Bower 1991), we take the 
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fiducial threshold mass of 10, 30 and 50 M* for the low, high 
and super-high accretion rate models, respectively. 

For a particular accretion history, the evolution of the 
morphological mix (fraction of E, SO and Sab galaxies) in 
clusters is calculated as follows: 

/ B (<) = - A 1 „, {A C um(tl) X f E {tl) + 



/so(t) 



^4cum (t) 
(Acum(t)- Acum(tl)) X /|° M }, 



{A C um(tl) X / S0 (*l) + 



(3) 



{A cum (t) - Aum(tl)) X /| 6ld + 
(-<4cum(£ — Ttrans) " J 4cum(^l — ' 7~trans)) 

x /sp 1 ^so}? 



(4) 



/Sab(t) = 7TT{(Aum(t)- J 4cum(t-r tr an S ))x/|p 1 ^ SO },(5) 

where t\ is the cosmic time at z = 0.54 (9.8 Gyr in our 
cosmology) the mean redshift of the four most distant clus- 
ters in the D97 sample, the epoch at which we normalise 
the fraction of ellipticals and SO's (/ E (ti) = 0.40 ± 0.08 and 
/so(ii) = 0.16 ± 0.04, see below). rtrans is the timescale for 
the morphological transformation from a spiral to SO after 
accretion, and /fp^so i s the fraction of field galaxies that 
are eventually turned into cluster SO's. Based on §2.2 this is 
taken as: 



f field 
JSp 



Id _ J Ji 

I Ji 



field 
Sab 



+ Js 



= 0.25 ±0.01 
= 0.55 ±0.02 



Model (c), 
otherwise. 



(6) 



In the following we also assume that any bulge-strong 
galaxies observed in the clusters are in the process of mor- 
phological transformation into SO's and we ignore merging 
of galaxies, since it is unlikely to be important for typical 
galaxies in these rich cluster environments at z g 0.5, where 
the relative velocities are high, Z 1000 km s -1 . 

To normalise the model for the evolution of the mor- 
phological mix within distant clusters we start from the 
morphologically-classified galaxy population seen within 
clusters at the highest redshifts covered by our analysis, 
z ~ 0.5. For this we use the average of fractions for the 
four most distant clusters in the Smail et al. (1997b) sample 
used by D97: CI 0412-65 (z = 0.51), CI 1601±42 (z = 0.54) 
C10016±16 (z = 0.55) and C10054-27 (z = 0.56). These 
are: E, 0.40 ± 0.08; SO, 0.16 ± 0.04; Sab, 0.32 ± 0.10; Scdm, 
0.07 ±0.06; and 0.05 ± 0.03 Irr/compact/unclassified. The 
errors show the scatter between the individual clusters. The 
Poisson errors are always small; about half of the errors 
shown. These numbers should represent the morphological 
mix for a 'typical' z ~ 0.5 cluster, since the fractions are sim- 
ply averaged across the four clusters rather than weighted 
by richness. 



4 RESULTS 

We plot in Figs. 1 and 2 the model S0/E and Sab/SO ratios as 
a function of redshift calculated from our model, normalised 
to the average cluster population at z = 0.54. Fig. 1 shows 
that the evolution of the S0/E ratio depends largely on the 
recent accretion rate, A TCCCnt , (compare Models (a),(b) and 



(f)) and/or what proportion of the accreted spiral popu- 
lation are turned into SO's (Models (a) and (c)), and that 
it is relatively insensitive to the transformation timescale 
,7"trans (Models (a),(d) and (e)). To match the strong evo- 
lution seen in the observations we must have both a high 
accretion rate, typically £ 50 per cent of the cluster popu- 
lation added since z ~ 0.5, and also transform the majority 
of the accreted spirals into SOs (not only bulge-strong Sab 
galaxies, but also bulge- weak Scdm's). We discuss the pos- 
sibilities for transforming bulge-weak spirals into SO's in §5. 
Based upon Fig. 1 we conclude that Models (a), (d), (e) 
and (f) are preferred. We note that if we applied the type- 
dependent magnitude cut for the accreted field galaxies to 
take into account the fading effect discussed above, the nec- 
essary accretion rate would need to be even higher since 
there would be fewer bright spirals available to be added 
from the field. This would strengthen the conclusions pre- 
sented here. This result confirms quantitatively the earlier 
discussion by Fabricant et al. (2000), who suggested that a 
high accretion rate of field galaxies was needed to explain 
the strong evolution in the morphological mix in clusters 
between z = 0.33 and the present-day. 

A further constraint on the models comes from the spi- 
ral fraction (Fig. 2). As our model assumes that most spiral 
galaxies observed in the cluster cores are undergoing a mor- 
phological transformation into SO's, the fraction of spiral 
galaxies in clusters can be used to estimate the timescale 
of this process (Poggianti et al. 1999). However, the spi- 
ral fraction also depends on the accretion rate and/or the 
proportion of the total spiral population which can trans- 
form into SO's. As the SO fraction depends upon the same 
parameters, we can reduce this dependence by normalising 
the Sab fraction using the SO fraction. The Sab/SO ratio 
thus depends primarily on the morphological transforma- 
tion timescale (Fig. 2) and can be used to constrain this 
timescale. 

We show in Fig. 2 the predicted Sab/SO ratios for our 
models and compare these to the observations. To illustrate 
the dependence on the transformation timescale we show 
three models: (d), (a) and (e) for which r tr ans = 3, 2, 1 Gyrs 
respectively. As Fig. 2 shows, if rtrans is too short (1 Gyr in 
Model (e)), the model struggles to retain enough accreted 
galaxies with spiral morphologies in the clusters before they 
are transformed into SO's. On the other hand, Model (d) 
(nrans = 3 Gyr) tends to over-produce the proportion of 
cluster spirals compared to the observations. Equally, adopt- 
ing a very high accretion rate model, Model (f), which is 
consistent with the S0/E ratio evolution (Fig. 1), and r tr ans 
of 1 Gyr tends to underestimate the Sab/SO ratios (Fig. 2). 
Taking into account the systematic uncertainty in the nor- 
malisation of the models at z ~ 0.5 using the fractions of E 
and SO galaxies in the four distant clusters (Fig. 1 and §3), 
our model suggests that the timescale for the morphological 
transformation should be in the railge Ttrans ^ 1-3 Gyr. 

Note that we have included field Scdm's into the cluster 
Sab category in the models, except for Model (c) (Eqs. (5) 
and (6)), since it is possible that the accreting Scdm's are 
transformed into Sab's either prior to, or immediately after, 
their entry into the cluster cores (§5). If we did not take into 
account this 'pre-transformation' from bulge-weak spirals to 
bulge-strong ones, the Sab/SO ratios in the models would go 
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(f) shi-all-1 
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Redshift (z) 
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Figure 1. The predicted evolution of the SO/E ratio in the 
cluster cores. The six lines show Model (a) [solid line]: high 
accretion, ttrans = 2 Gyr, Model (b) [dotted line]: low accre- 
tion, ttrans = 2 Gyr, Model (c) [dashed line]: high accretion, 
ttrans = 2 Gyr, Model (d) [long dashed line] : high accretion, 
ttrans = 3 Gyr, Model (e) [dot-dashed line]: high accretion, 
ttrans = 1 Gyr, and Model (f) [long dot-dashed line]: super- high 
accretion, ttrans = 1 Gyr, respectively. The morphological trans- 
formation of Sab+Scdm — ► SO is assumed, except for Model (c) 
where Sab — > SO is assumed. Observed ratios are taken from D97 
(circles), C98 (triangles) and D80 (squares). The average ratio in 
the D80 sample and scatter (lcr) are shown by the filled square 
and the error-bar respectively. Note, the systematic uncertainty 
in the normalisation of the model based on the cluster popula- 
tions at 2 = 0.54 is ±0.15 in the S0/E ratio and this produces a 
systematic uncertainty in the S0/E ratios at z ~ of 15-20%. 



Figure 2. The evolution of the Sab/SO ratio in our model clus- 
ters. The four lines show Model (d) [long dashed line]: ttrans = 
3 Gyr, Model (a) [solid line]: f tra n S = 2 Gyr, Model (e) [dot- 
dashed line]: ttrans = 1 Gyr, and Model (f) [long dot-dashed line]: 
ttrans = 1 Gyr, respectively. The high accretion rate is assumed, 
except for Model (f) where the super-high accretion rate is as- 
sumed. The data points are taken from D97 (circles) and C98 
(triangles) and D80 (squares). The average and scatter (lcr) are 
shown by the filled square and the error-bar, respectively. Note 
that the uncertainty in the model normalisation at z = 0.54 using 
the E and SO fractions produces a systematic uncertainty in the 
Sab/SO ratios of 30-35% at 2 = 0.54, 10-15% at z ~ 0.4, and 
~4% at z ~ 0. 



down by ~ 50 per cent. In this case, only Model (d) would 
reproduce the observed trend, suggesting r tmns ~ 3 Gyr. 

The relatively long timescale required for the morpho- 
logical transformation in the model is consistent with a num- 
ber of recent results on the properties of galaxies in distant 
clusters. Based upon the paucity of post-starburst signa- 
tures in the spectra of the luminous SO population in dis- 
tant clusters and the prevalence of red, passive galaxies with 
late-type morphologies, Poggianti et al. (1999) concluded 
that the morphological transformation of SO's had to oc- 
cur on a longer timescale than the stellar evolution lifetimes 
(Z 1 Gyrs) of any A-stars formed in the most recent star 
formation event in these galaxies (either while they were in 
the field, or any activity which was induced on entering the 
cluster). Comparisons between the properties of the lumi- 
nous SO and elliptical galaxies in distant clusters have been 
used to place limits on the relative evolution of the two pop- 
ulations: work on three clusters at z = 0.31 - using spectral 
(Jones, Smail & Couch 2000) and photometric (Couch et al. 
in preparation) analysis concluded that the SO and elliptical 
populations are remarkably homogeneous, and similar con- 
clusions were reached using precise photometric analysis of 



the colours of SO/ellipticals in the cluster cores at z = 0.18 
(Smail et al. 2001), z = 0.33 (van Dokkum et al. 1998) and 
z ~ 0.5 (Ellis et al. 1997). All of these studies suggest that 
the luminous SO galaxies in the cluster cores ceased their 
star formation at least 2 Gyr prior to the observed epoch 
(but see the discussion in Smail et al. 2001 for less lumi- 
nous SO's). Observations of the evolution of the colour dis- 
tribution of cluster galaxies (Kodama & Bower 2001) and 
the radial gradients in star formation, measured using [On], 
(Balogh, Navarro & Morris 2000) both suggest that the star 
formation in a galaxy declines relatively slowly after it has 
been accreted by a cluster, ~ 1 Gyr. 

Taken together these results would indicate that an ac- 
creted star-forming spiral galaxy transforms into a passive 
SO within <; 2 Gyrs of entering the cluster environment. Our 
analysis confirms this suggestions. This is encouraging as 
our analysis relies purely on modelling the evolution of the 
morphological mix within the clusters, whereas the previous 
results are based on the stellar populations within individual 
cluster galaxies and so these two approaches are completely 
independent. 
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5 DISCUSSION AND CONCLUSIONS 

We illustrate a simple phenomenological model to attempt 
to reproduce the broad properties of the morphological evo- 
lution of cluster populations across ~ 7Gyrs from z ~ 0.5 
to the present-day. 

In the framework of this model, we find that the ob- 
served evolution of the SO/E ratio requires high rates of 
accretion onto the clusters and also that a significant frac- 
tion of all accreted spiral galaxies, both bulge-strong Sab 
and bulge-weak Scdm, are transformed into SO's. 

The need to transform not only Sab, but also weaker- 
bulged galaxies into SO's is a concern as there is no obvious 
and effective mechanism for transforming an Scdm into an 
SO. The bulges of late-type spirals are just not luminous 
enough to evolve into the bright bulges of SO galaxies sim- 
ply as a result of stopping the star formation in the disk or 
even by 'harassment' in the cluster environment (Moore et 
al. 1996). One possibility would be that the accreted field 
populations are not the same as 'pure' field population we as- 
sumed from the MDS. They may be 'pre-processed' in group 
or sub-cluster units, or intrinsically biased towards earlier 
types. Most of the bulge-weak Scdm's may have been al- 
ready turned to bulge-strong galaxies by merging or through 
bar- formation due to strong galaxy harassment/interaction, 
before entering the clusters. This is likely, given that the 
outskirts of a rich cluster, where accreting spirals are lo- 
cated, also originate from high density regions in the early 
Universe (e.g. Moore et al. 1999b). Therefore, the morpho- 
logical mix in the field local to the cluster may be skewed 
towards earlier Hubble types than that in the 'pure' field 
through galaxy-galaxy interactions in these denser regions 
(Moore et al. 1999a, 1999b; Springel 1999; Diaferio et al. 
2001). 

An alternative solution is to adopt a high density Uni- 
verse. With Q. — 1 (qo = 0.5) for example the accretion rate 
is higher than in the current models A TCCCnt (z = 0.5) ^> 0.7 
(Bower 1991). And if this is really the case, the large S0/E 
ratio at z ~ can be reached without the need to trans- 
form Scdm galaxies into SO's. Furthermore, the SO/Sab ratio 
would also be consistent even with r tra ns = 1 Gyr. However, 
recent observations suggest slower cluster evolution than ex- 
pected in an 0, o = 1 cosmology, with little evolution seen in 
X-ray luminosity function of clusters out to z ~ 0.3 (Ebcling 
et al. 1997) and beyond (Fairley et al. 2000) , indicating the 
recent accretion rate onto rich clusters cannot be that high. 

Finally, we note that there is the possibility of signif- 
icant scatter in the morphological mix in clusters even at 
the same epoch (Figs. 1 and 2). However, with the present 
observations this scatter is consistent with the observational 
errors and so we have chosen to model the evolution of the 
morphological mix in an 'average' cluster. If real differences 
do exist between the morphological distributions in similar 
clusters at a fixed redshift this would provide a powerful tool 
for investigating the nature of galaxy morphology and the 
mechanisms responsible for its transformation. 

In summary, we have simulated the evolution of the 
morphological mix in clusters with the aim of applying a 
phenomenological test for the hypothesis of the morphologi- 
cal transformation. We therefore use a simple model in which 
the accreted field spirals are transformed into SO's in the 
cluster environment to attempt to match the observed evo- 



lution of the distribution of morphologies for cluster galaxies 
out to z ~ 0.5. To reproduce the rapid increase in the S0/E 
ratio between distant clusters and their present-day counter- 
parts, it is found that at least nearly half of the present-day 
cluster population must have accreted since z ~ 0.5, con- 
sistent with the earlier estimate by Fabricant et al. (2000). 
Moreover, we found that for this hypothesis to be correct, 
not only the bulge-strong field spirals, but also relatively 
bulge- weak field spirals (Scdm) must have eventually trans- 
formed into cluster SO's. This requirement may be difficult to 
achieve unless the field population is pre-processed to pro- 
duce stronger-bulged galaxies prior to their arrival in the 
cluster core. Finally, in order to explain the relatively high 
fraction of cluster spirals we find that the process responsi- 
ble for the morphological transformation from spiral to SOs 
upon entry to clusters appears to be relatively slow (~1- 
3 Gyr). 
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